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We study the charmful three-body baryonic B decays with _D'*' or J/^ in the final state. 
• We explain the measured rates of — > npD*^ , & — > ppD^*^^ , and _B~ — > KpJ In 

^1^' particular, the branching fraction of AAD" predicted to be of order 2.3 X 10~® is 

' \ in accordance with the Belle measurement, (1.05tQ44 =•= O-l^) ^ < 2-6 x 10~^. 

I) . 

. 1. Introduction 

The experimental measurements present that the charmful and charmless three- 
body baryonic B decays have the same features. First, the threshold effect is ob- 
served in 5° ppD^*^^,^ where a curve peaks near the threshold area in the 
CO . dibaryon invariant mass spectrum. This phenomenon found in all charmless cases 

' can be understood in terms of a simple short-distance picture.^ Particularly, this pic- 



O 



X 



ture has been used to realize why charmless three-body decays^ ^ have rates larger 
^ ! than charmless two-body decays;i°'" that is, T{B BB'M) > T{B BB'). One 

' energetic qq pair must be emitted back to back by a hard giuon in order to produce 

■ a baryon and an antibaryon in the two-body decay. This hard gluon is highly off 

^ I mass shell and hence the two-body decay amplitude is suppressed by order of as/q^. 

In the three-body baryonic B decays, a possible configuration is that the BB' pair 
is emitted coUinearly against the meson. The quark and antiquark pair emitted 
I from a gluon is moving nearly in the same direction. Since this giuon is close to 

its mass shell, the corresponding configuration is not subject to the short-distance 
suppression. This implies that the dibaryon pair tends to have a small invariant 
mass.^^ 

Second, the Dalitz plot of i? — > ppD^*^^ with asymmetric distributions signals 
a nonzero angular distribution asymmetry as measured in B^ — )■ ppn^ , B^ 
A+p7r+, B~ Ap7, B~ ppK~, and B~ Apir^ .^^^-^"^ Nonetheless, while 
the above-mentioned short-distance picture predicts a correct angular distribution 
pattern for the decays PpT^ , B^ A+p7r+, and B^ ^Fl, it fails to 

explain the angular correlation observed in B^ —>■ ppK^ and B^ Kpn^ . The 
intuitive argument that the K~ in the pp rest frame is expected to emerge parallel 
to p is not borne out by experiment.^"* Likewise, the naive argument that the pion 
has no preference for its correlation with the A or the p in the decay B~ ApTT~ is 
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ruled out by the new Belle experiment^ in which a strong correlation between the 
A and the pion is seen. 

Therefore, the study of the charmful baryonic B decays B BB'Afc may 
help improve our understanding of the underlying mechanism for the threshold 
enhancement and the angular distribution in three-body decays. We shall therefore 
focus on B ^ BB'Mc with = D^*^ or J/* to see if we can explain the branching 
fractions. 

2. Experimental Data 

Table 1. Branching ratios of iJ — ► BB'A/c decays in units of 10"* for Mc = Z)'*' and 
10^'' for Mc = J/^. 



Decay CLEO^^ BaBari'i^ Belle"-!'' 



50 


npD*+ 


14.5t|;^ ± 2.7 




B° 


ppD° 


1.13 ±0.06 ±0.08 


1.18 ±0.15 ±0.16 


50 


ppD*° 


1.01 ±0.10 ±0.09 


1.20lo Jg ±0.21 


B- 


— > ppD~ 




< 0.15 


B- 


ppD*~ 




< 0.15 


B° 


-^ApD+ 




0.29 ±0.07 ±0.05 ±0.04 


§0 


^ AADO 




{^■05^0,11 ± 0.14 < 2.6) X 10-1 


B- 


^ ApJ/* 


11 6+'''-^+'*-^ 


11.6 ± 2.8+2;3 


B- 


TPpJ/-^ 




< 11 


B° 


ppj/'^ 


< 1.9 


< 0.83 



The branching fractions of the charmful baryonic B BB'Afc decays are 
summarized in Table [H where — > npD*^ was first observed by CLEO in 
2001.^^ Note that the decays ppD° and 5° ppD*" have similar re- 

sults in rates. The nonobservation of B^ — > ppD'-*^^ is due to the fact that 
it proceeds via b —^ ucd at the quark level. This leads to a suppression of 
l^u&Kd/^c&K,*;;!^ — 10^'' compared to its neutral partner. Likewise, it is expected 
that B{B° ppj/^) = \Vcd/Vcs\'^ Br{B- ApJ/^) ~ 10"^, consistent with the 
experimental upper bound for this decay mode. As for the rate difference between 
B'-^ ppD^ and B'^ AAD'^, it has to do with the baryonic form factors, which 
we are going to elaborate on later. 

3. Formalism 

To have the amplitudes, we shall adopt the generalized factorization approach, ^"'^^ 
which has been well applied to the study of three-body baryonic B decays. ^"^'^^"^^ 
Under the factorization approximation, the decay amplitudes can be classified into 
three different categories: the current-type (class-I), the transition-type (class-II), 
and the hybrid-type (class-Ill) amplitudes. The amplitude of 5° npD*~^ is the 
class-I current-type which proceed via a color-allowed, external T4^-emission diagram 
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(a) (b) 
Fig. 1. Two types of the B —* BB'A/c decay process: (a) current type and (b) transition type. 



as depicted in Fig. [TJa) , which is given by 
AciB" npD*+ 



^VcbV:,af {np\idu)v-A\0){D*+\{cb)v^A\B°) , 



V2 



(1) 



with af to be specified later. The amphtudes of -> ppD^*^^ , AAD°, 
B^ Ap(S°p) J/4', and ppj /'^ are the class-II transition- type via the color- 

suppressed internal W emission diagram [Fig.[TJb)]. The amplitudes read 



Ar{B ^ BB'M,) = ^V,,V; 



Mr 



{M,\{cq)v-A\0){BB'\iq'b)v-A\B) , (2) 



with a2 " to be given later, where B BB'A/c could be B° ppD'^*'>°, AAD° for 
qq' = ud, B^ ~f ppj/^ for qq' = cd, B^ ApJ/^, YPpJ/^ for qq' = cs. 
For the dibaryon creation in Eq. ([T]), we write 



(BB'|gi7^g2|0) 
(BB'|gi7^75(Z2|0) 



hA{t) 



(3) 



where u{v) is the (anti-)baryon spinor, and i^i.2, Qa, hA are timelike baryonic form 
factors. Note that there are two additional form factors in the form of uq^^v and 
ua^v<lul5V. However, since we assume SU(3) flavor symmetry, we can neglect these 
two form factors as they vanish for conserved currents. The asymptotic behavior 
of form factors is governed by the pQCD counting rules. '^^'^^ In the large t limit, 
the momentum dependence of the form factors Fi (t) and {t) behaves as 1 /t'^ as 
there are two hard giuon exchanges between the valence quarks. More precisely, in 
the t oo limit 



9A{t) = ^ 



(4) 



where 7 = 2-1-4/(3/3) 2.148 with /3 being the QCD (3 function and Aq = 0.3 GeV. 
In the asymptotic t 00 limit, both F2{t) and hA{t) have an extra l/t dependence 
relative to Fi and qa owing to a mass insertion at the quark line.'^"''*'^^'*^ However, the 



form factor Ha is related to gA by the relation of Ha — — 5a(™b 



/)^/t, through 



the equation of motion. Hence, in ensuing numerical analysis we will keep /lyi(t) and 
neglect F2{t). Under the SU(3) flavor and SU(2) spin symmetries, the parameters 
Cpi and Cg^ appearing in — > BB' transitions arc no longer independent but are 
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related to each other through the two reduced parameters C|| and Cjj. Then we 
have 

Cf, = - ICji , Cg, - + ICj^ , for -> np, (5) 

As for the three-body transition B BB', its most general expression reads 

{BB'\Vf,\B) = m[5i7^ + yaicr^i/P" +33^^ + 54(^6' +Pb)p + gsiPw -Pb)m]75W, 

(BB'I^^JB) = m[/i7p + /2icr^./p'' + /sP^ + .U{Pb' +Pb)p + /sbg' -Pb)p]w,(6) 

with V{A)fi = q'jfj,{'j5)b and p — pb — P'b — Pu'- Since three gluons are needed 
to induce the B — > BB' transition, two for producing the baryon pair and one for 
kicking the spectator quark in the B meson, the pQCD counting rules imply that 
to the leading order 

Mt) = ^ , mit) = (7) 

Just as the previous case for vacuum to the dibaryon transition, under the SU(3) 
flavor and SU(2) spin symmetries, the parameters Dg-^, ^/u Dg. and Df. can be 
expressed in terms of the reduced parameters D^, Dj^ and Z3^|. Then we have 

DPP(DPP) = T iDj^, DPP = -Df ^ -^D'l^ , for B° pp, 

KHd^^) - 5^11 + h^Jl. Df^ = Df = 0, for 50 ^ AA, 

^gf Kf ) = -^i^lh D^! = -I^U=-^2D\v ioi-B-^Ap, (8) 
where the upper (lower) sign is for Dg^ {Df^) with i — 2, 3,..., 5. For the meson 
parts, the values for the decay constants of and J/^E* can be found in Refs. 

43,44. The definition and parametrization of B to Z?'-*-' transition form factors can 
be found in Refs. 45,46. 

4. Numerical Analysis 

We need to specify various input parameters for a numerical analysis. For the CKM 
matrix elements, we use the values of the Wolfenstein parameters in Ref. 47. Here, we 
extract af*, af**\ and 03^* from npD*+ , B° ppD*", and B- ApJ/^, 

respectively, which are given by 

af = 1.23 + 0.19 , af'*' = 0.33 + 0.04 , as^* = 0.17 + 0.03 . (9) 

We note that factorization works if the parameters ai and 02 are universal; namely, 
they are channel by channel independent. Since af^" lie in the ranges of ai ~ 0{1) 
and a2 ^ 0(0.2 — 0.3),^^ which is suggested by two-body mesonic B decays, we 
shall assume the validity of factorization in charmful baryonic B decays. 

For the parameters C|| and Cjj in Eqs. ([31 U]), we use the data of e+e^ 
pp, nn^^''^^ to determine their magnitudes and the decay rate of B^ npD*^ to 
fix their relative sign 

(q,, = (67.9 + 1.4, -216.9 + 23.5) GeV . (10) 
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Table 2. Branching ratios of _B — > BB'Mc decays, where the first and 
second errors come from the uncertainties of ai (02) in Eq. Q and 
baryonic form factors in Eqs. IIIOI llll l. respectively. 



B X 10-'* 



5° -> npD*+ 
B° ppD° 



Our work 



14.4 ±4.8 ±3.2 
1.1 ±0.3 ±0.2 
1.0 ±0.3 ±0.4 



B X IQ-'^ 



B° -> AAD" 

B- T,°pJ/'i> 
B° ppj/^! 



Our work 



2.3 ±0.6 ±0.6 
11.6 ±4.5 ±2.9 
0.13 ±0.05 ±0.02 
1.2 ±0.5 ±0.2 



As for the parameters D\\ and Uyj- in Eqs. ([61 [7]), we employ the observed rates 
of 5° ppD^, B- ppK*-, B° ppK*°, and B~ pp^^-ifi^Q^n, 50-52 
conjunction with the measured angular distribution of the last decay mode to obtain 

(Dll, D^) = (67.7 ± 16.3, -280.0 ±35.9) GeV^ 
(DlDl,Dl,Dl) = 

(-187.3 ± 26.6, -840.1 ± 132.1, -10.1 ± 10.8, -157.0 ± 27.1) GeV' . (11) 

To calculate the decay rates, we use the equation in Ref. 43 for three-body de- 
cays. We present the numerical results for the branching ratios in Table [21 and the 
dibaryon invariant mass spectrum for i?° — > ppD'^*^^ and B^ — > ApJ/^ in Fig. [51 
We note that the first and second errors in Table [H come from the uncertainties of 
ai (02) in Eq. ^ and baryonic form factors in Eqs. pTD fTT|) . respectively. 



tS 2.5 

\2.0 

I 0.5 
0.0 



2.5 3.0 



~3.0 
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400- 

>: 

^300 

55 
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Fig. 2. Dibaryon invariant mass distributions for B*^ — > ppD'^, — ► ppD*^ and B — > A.pJ 
respectively. Experimental data are taken from Refs. 1,19. 



5. Discussion and Conclusion 

As seen in Table [21 our prediction B{B^ Yj^pJ/^>) ~ 1.3 x 10"^ is consistent 
with the Belle limit, 1.1 x 10"^, and B{B° PpJ/^) = 1-2 x 10"'' is in accordance 
with the BaBar limit but slightly higher than the upper bound set by Belle. As for 
the threshold peaking effect, while it manifests in the decay i?" ppD'^*^^ the data 
clearly do not show the threshold behavior in B~ ApJ/"^ (see Fig. [2]). This can be 
understood as follows. In the latter decay, the invariant mass mAp ranges from 2.05 
to 2.18 GeV, which is very narrow compared to the nipp range in the B° ppD^*'^° 
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decay. Consequently, the invariant mass distribution of dT/dm/^p is governed by 
the shape of the phase space due to the relative flat dependence within the 
small allowed mAp region. B{B'^ ~> AAI?") predicted to be of order 2.3 x 10^^ is 
consistent with the measured Belle data, (imtoil ± 0.14) x 10"^ < 2.6 x 10"^ 
B{B^ AAD'^) < B{B'^ ppD^) can be understood by the baryonic form factors 
in Eqs. ((8l fTTj) . which are evaluated to be 

[Df^, Df^) = (172.8, -105.1) GeV^ , {Df, Df^) = (207.9, -162.7) GeV^ , 

{D^^^, Df) - , {Dll Df^) ^ . (i-2, 3, 5) (12) 

In sum, within the framework of the generalized factorization approach, we have 
explained the measured charmful three-body baryonic B decays with or J/^P in 
the final state. The measured B{B° AAD^) < S(i3" -> ppD°) can be understood. 
This is due to the constraint in the baryonic transition form factors for B^ AA 
in the approach of the pQCD counting rules. 
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